Cooper channel and the singularities in the thermodynamics of a Fermi liquid 
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We analyze how the logarithmic renormalizations in the Cooper channel affect the non-analytic 
temperature dependence of the specific heat coefficient -f{T) — 7(0) = A(T)T in a 2D Fermi liquid. 
We show that A(T) is expressed exactly in terms of the fully renormalized backscattering amplitude 
which includes the renormalization in the Cooper channel. In contrast to the ID case, both charge 
and spin components of the backscattering amplitudes are subject to this renormalization. We show 
that the logarithmic renormalization of the charge amplitude vanishes for a flat Fermi surface, when 
the system becomes effectively one-dimensional. 

PACS numbers: 



Non-analytic behavior of the thermodynamic param- 
eters of a Fermi liquid has a been of interest since the 
early days of Fermi-liquid theory, when it was found that 
for D = 3 the specific heat coefficient 7 (T) = C(T) /T 
has a non-analytic temperature dependence Sj(T) = 
7(T) - 7(0) oc T 2 InT [if. The issue of non-analyticities 
has been revived in recent years, following the work by 
Belitz, Kirkpatrick, and Vojta [2|, who found that the 
same type of logarithmic behavior holds for the spin sus- 
ceptibility at a finite momentum q in three dimensions 
(3D) 5 Xs (q)(xq 2 lnq. 

Non-analyticities are stronger in two than in three di- 
mensions: (J7 (T, H) and Sxs{T\ QaJQ, are linear functions 
of their respective variables [EM IE IE 0, IE IE EE H 
H EE El El- In Refs. 0, El, it was argued that 
the non-analytic part £7 (T) comes from one-dimensional 



scattering processes embedded in a 2D phase space [14 1 
and that Sj(T) is expressed via the charge and spin com- 
ponents of the exact scattering amplitude with zero total 
momentum f(9 — tt) ("backscattering amplitude") 



3C(3) 
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[f c (n)+3f 2 (n)]T. (1) 



Validity of Eq.([T|) was verified perturbatively in Ref. 11 1 
by calculating 5~/(T) and / CjS (7r) independently to third 
and second order in the interaction, respectively, and 
checking that Eq.(TT]) holds. That analysis, however, 
was incomplete - it included the renormalizations in the 
particle-hole channel but neglected the renormalizations 
in the Cooper channel 



Rcfs. 



Ill I13 ] conjectured- 



without proof- that Eq. |T]) still holds if the Cooper renor- 
malizations are included into f c , s (ir). 

The Cooper renormalization of the backscattering am- 
plitudes is an essential ingredient of the theory, partic- 
ularly in the limit T — > 0. The argument is that the 
backscattering amplitude 

is expressed in terms of two full vertices with zero to- 
tal momentum and cither zero or 2kp momentum trans- 



fer r(k,-k;k,-k) 
r(k,-k): 



T(k,k) and r(k,-k;-k,k) 



/cW = 



r(k,k)--r(k,-k) 



,/ s (tt) = -— r(k,-k) 

(2) 

Equivalently, each of these two vertices can be ex- 
pressed via the fully renormalized Cooper vertex with 
either zero or 2kp momentum transfer: T(k, k) = 
r c (0), r(k,-k) = T c (2k F ). Therefore, these vertices 
can be expressed via the partial components J„ of the 
irreducible interaction in the Cooper channel as |20] 



r c (o),r c (2fc F ) = ^(±i 
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At low temperatures, all terms in the sum scale as 
(IuEf/T)- 1 , so that T(k, k) andT(fc, -k) should, in gen- 
eral, be reduced by the same logarithmic factor. Whether 
this logarithmic renormalization affects the backscatter- 
ing amplitudes and Sj(T) is a more subtle issue which is 
to be addressed by a direct calculation. 

The ID case serves as a good example here. Both 
r(fc, k) and T(k, —k) are logarithmically renormalized in 
ID, yet, these renormalizations only affect the spin chan- 
nel, but cancel out in the charge channel IE 12, 18, IE]- 
As a result, the specific heat remains linear in T at the 
lowest T. This agrees with ID bosonization according to 
which charge excitations are described by a free Gaussian 
theory, whereas the spin channel for the case a repulsive 
interaction contains a marginally irrelevant perturbation 
that causes a logarithmic flow of the spin amplitude. The 
issue that we address in this paper is whether this situa- 
tion occurs only in ID or in higher dimensions as well. 

The interplay between the logarithmic renormaliza- 
tion of the interaction and the behavior of specific 
heat coefficient in 2D has recently been considered by 
Aleiner and Efetov (AE) [IE] (for a subsequent analy- 
sis see [E 111)- AE invented an elegant supersymmet- 
ric method to treat the problem in arbitrary D by inte- 
grating out fermions and expressing the low-energy ac- 



tion solely in terms of the low-energy collective bosonic 
modes. They found that the spin contribution to 5j(T) 
is affected by the Cooper renormalization and behaves 
as T(ln|lnT|/lnT) 2 in the limit of T -> 0. They 
treated the charge component in the eikonal approxima- 
tion which neglects the curvature of the Fermi surface, 
and found that, in this approximation, the charge com- 
ponent remains unrenormalized. 

In this communication, we re-consider this issue. We 
evaluate 5"f(T) explicitly to third order in the interaction 
U(q), including Cooper renormalizations, and also eval- 
uate the spin and charge components of the backscatter- 
ing amplitude to second order in U(q). We find that, in 
contrast to the ID case, both / s (7r) and / c (tt) undergo 
logarithmic renormalization in 2D. We also find that Eq. 
((T|) still holds when U(q) is re-expressed in terms of / s (7r) 
and /c(tt)- This result agrees with the conjecture made 
in Refs. [ll|, [l2j]- For a short-range interaction, we find 
that the spin contribution to 5~f(T) scales as (1/lnT) 2 , 
while the charge contribution scales as (In | lnT|/ InT) 2 . 
As a result, if the system remains in the normal state 
down to very low T, Sj(T) would be dominated by charge 
fluctuations and behaved as 5~f(T) oc T(ln | lnT|/ InT) 2 . 
The normal state behavior in 2D, however, exists only 
down to the temperature T p of the Kohn-Luttinger in- 
stability towards p— wave pairing 2^, 23[. We analyzed 
5 1 (T) = A(T)T near T p and found that A(T) diverges 
as(T-T p )- 2 . 

As a simple check, we verified that the same procedure 
that we used for calculating the scattering amplitudes 
in 2D reproduces the known result in ID, namely, the 
cancellation of the lo gar ithmic renormalizations in the 
charge channel TfJ This comparison helps to 

see where exactly the cases of D > 1 and D = 1 dif- 
fer: for D > 1, only the Cooper channel is logarithmic; 
hence there is no cancellation between the Cooper and 
particle-hole channels, and f c (n) is renormalized along 
with f s (it) . For D = 1, both particle- hole and particle- 
particle renormalizations are logarithmic, and the two 
cancel each other in the charge channel. To further 
emphasize this point, we also considered the 2D case 
with a non-circular Fermi surface and demonstrated that 
when the Fermi surface becomes fiat, the particle-hole 
renormalization becomes logarithmic and the Cooper and 
particle-hole renormalizations again cancel each other in 
the charge channel. 

We believe that it is the effect of the curvature that is 
responsible for the difference between our result and that 
of AE. 

We begin with the 2D case and a circular Fermi surface. 
To second order in U(q), the diagrams for the thermody- 
namic potential 5 contain two fermionic bubbles Tl(q, ft). 
To this order, the non-analyticity comes from the square 
of the dynamic part of n((/, f2) which contains the term 
fl 2 /q 2 at VFq 3> | SI | 0, [U, EH- Because the momentum 



integral is logarithmic in 2D, the second-order thermo- 
a b / N Y c 







FIG. 1: Diagrams for the thermodynamic potential to third 
order in the interaction. 



dynamic potential 



d 2 qU 2 ~ [/ 2 T^S! 2 ln|S!| 

n 



U 2 T 3 (4) 



contains a universal T 3 term, which gives rise to a 0(T) 
term in 5j. It was shown in Refs. 0, [HI, that 
the momenta carried by fermions in the two bubbles, 
ki,k.2, k3, and k4, are correlated in such a way that 
ki « k.2 ~ — k,3 w — k4. These four fermions can then be 
re-arranged either into a convolution of two particle-hole 
bubbles H p h(fl,q) or two particle-particle bubbles @, 
H P p(tt,q); the term Q 2 /q 2 is produced regardless of the 
choice. 

Relevant third-order diagrams for S are shown in Fig[T] 
They contain either three particle-hole or three particle- 
particle bubbles. It can be shown that the non-analytic 
terms in Sj(T) come from the terms containing the prod- 
ucts of two dynamic and one static parts of these bub- 
bles. Two dynamic parts of the bubbles produce fl 2 /q 2 
term, which is the source of non-analyticity, whereas the 
third bubble renormalizes static backscattering vertex. 
If the third bubble is a particle-hole one, the renormal- 
ized vertex is a constant. The momentum integration 
then yields fl 2 In | fl \ , and subsequent frequency summa- 
tion gives rise to the T 3 -term in H. In the third bubble is 
a particle-particle one, the renormalized vertex contains 
an additional factor of In q. This changes the result of the 
momentum integration to fl 2 In 2 [note an extra In |f2| 
as compared to Eq. 0$] and gives rise to a T In T-term in 
Sj(T). 

The contributions from diagrams (a-d) have already 
been presented in fill ], diagrams (e) and (f) were not 
considered there. For completeness, we present the re- 
sults for all third-order diagrams. We have [2l[ 
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S 3q = - (u o («0Utf_e) + 2u o u 7I (u ) + Un({ug)) + 2u u 7I ((ug)}) K, ^ 3b = (4ul(ug) + 2u§u,r) if 

+ [4t4((«s» + 2t4«o((l))] K, E 3c = -4 [«§ + t4((l»] if, S 3d = 2u 7r (u e u, r _ e )2<' + 2u ((i4»#, 

Ep Ep 
S 3e = -(2u (ue) + 2u 1T (u e u^-e))K\n— , S 3/ = (u^ugu^-g) + u v (uj))K In ■ 



T 



r 



(5) 



wherein = f(3)T 3 /7ruf. and u# = (m/2Tr)U(2kFSm0/2), Quantities denoted as (gg) and {(ge)) are (gg = ug,Ug, 
such that = (m/27r)£/(0) and u w = (m 12^)11 (2k p). and 1): 



2tt /-d 2 /^ ..,..^2 rde 2^ fd 2 ^ r<j6 



w d6> 9, , 

cos - In cot , 
2 V 4 



r 



For a circular Fermi surface, ((1)) = 1. w 2 — uqU^)K [ll[ we obtain a complete result for Sj(T) 

Combining S 3 with the second-order result S2 = (ujj + to third order in the interaction 



J 



Sj(T) = - Uq + u 2 n - uqUk - 4ul + 2uqU„ + (4u^ - 2u u v )(ug) - (tt - 2u 7T )(ugu 7T -g) + (-4?4 + 2it *4)((l)) 



+ (4m 2 - 2u u 7T )((u g )) - (u v - 2w )((u 2 ,}} - (2tt (ttg) + 2u n (ugu 7T -g) - u (ugu v -e) - u w (ttj)) In 



"t" 



6C(3)T (7) 



r 



Next, we compute independently the static interaction 
vertices r(k,k) and r(k, — k) (r fc in the Fermi liquid no- 
tations) to second order in U(q), including the renormal- 
izations in the particle-hole and particle-particle chan- 
nels. Collecting both contributions, we find 



£r(k,k) = ^0 + <(^1» - <^|)i n ^; 



fr ( k,-k) 

Z7T 



Ep 

(ugllTt-g) ( 1 - In — 



(8) 
(9) 



where 



ii = uq (1 - 2u + 2{ug)) , 

u„ = u„ (l-2u w ((l»+2((« fl »). 

The charge and spin components / c ,s(t) are obtained us- 
ing Eq. |J2J|. Evaluating the expression / 2 (7r) + 3/ 2 (7r) 
to third order in U(q) and substituting the result into 
Eq.([T]), we find that Eq. ([7]) is fully reproduced. There- 
fore, at least to third order, the prefactor of the T term in 
5j(T) is expressed via exact scattering amplitudes, which 
include the renormalizations in the Cooper channel. 

At low T, the particle-particle renormalizations are 
more relevant that those in the particle-hole channel, as 



the former contain a large factor of InEp/T. Keeping 
only the particle-particle renormalization, we obtain 



f c (n) = 2u -u n + ((ugiiTr-g) - 2{u 2 g )) In ^ 



f s (n) = - ( Utt - {ugu-x-g) In 



T 



(10) 



We see that both f s (n) and / c (tt) contain logarithmic 
corrections. For ug =const= u, Eqs. (jTDJ) reduce to 



/ c (7r) = —fs(n) = u I 1 - nln 



T 



1 + u In ^f- 



(11) 



A better estimate is obtained if we assume, following AE, 
a simple model form for the irreducible interaction in the 
particle-particle channel: J 

J c (q) = T(k, -k; k + q,-k-q) = aw/{q 2 + a 2 ), 

where q = 2fcp-sin#/2 and w has the units of velocity. 
The partial components of J c (q) are 



c 



we 0n „, / a / a 2 

— r ^ 21n fe + v + ^ 



\/a 2 + 4fc? 



(12) 
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In the limit of T — > 0, the sums over n in Eq. ([3]) for 
r c '(0) and r c (2fcir) are dominated by large n; replacing 
summation over n by integration, we obtain 



r , x 2tt InL r , s 2?r 1 - e"' 3 



(13) 



where L = \sxEp/T. In this limit, r c (0) is larger than 
r c (2M by InL, hence / c (tt) « 21nL/(/3L) > / s (tt), 
and the full result for the specific heat coefficient becomes 



corrections: 



««*» 



Lp 

u In— , 



(16) 

where now u# = U(2kp sin9 /2) / (2irvp). Combining the 
logarithms in the particle-particle and particle-hole chan- 
nel and neglecting non-logarithmic second-order terms, 
we reproduce the well-known results for the charge- and 
spin scattering amplitudes in ID [H, E3, El : 



&y(T) = - 



3C(3) 
2tt (v* F f 



2 In In %f 

/31n " 



E F 
T 



T. 



(14) 



Note that the prefactor depends on the functional form of 
r (q) but not on the magnitude of the interaction. Eq. 
(fl~3|) is valid only at such low temperatures that InL ^> 
1. For a more realistic case of L > 1, both r c '(0) and 
T c (2k F ) are of order 1/L, and 5-y(T) oc T/ (lnL F /T) 2 . 

A more fundamental reason why the ultra-low T 
regime is unaccessible is the Kohn-Luttinger effect: the 
superconducting instability for a nominally repulsive in- 
teraction [22l |. It is known that irreducible J c \q) is non- 
analytic near q = 2kp due to screening of the origi- 
nal pairing interaction by the particle- hole excitations. 
Screening produces a long-range component of J c (q) . In 
2D, the Kohn-Luttinger effect is a bit tricky as one needs 
to include vertex corrections to the polarization bub- 
ble to obtain an oscillating long-range component of the 
pairing interaction sin(2fc F r)/r 2 between particles at the 
Fermi surface [23| . Because of oscillations and 1 jr 2 be- 
havior, the partial harmonics J% acquire negative parts 
that fall off with n algebraically rather then exponen- 
tially: Jn' KL ~ -a/n 2 , a > (at small U, a oc U 3 ). 
As a result, J„ become negative for n > n c , which 
implies superconductivity. For a moderately strong in- 
teraction [mU(q) > 1], Jn becomes negative already 
for n — 1 (|Jf I is the largest), and the system be- 
comes unstable towards p— wave pairing at T p defined 
by (ma /2tt) In E F /T P = 1. Near T p , the n = 1 term 
dominates the sums in Eq. (j3J), both r c (0) and T c (2kp) 
diverge as a/(l — (ma/2n)L) oc 1/(T — T p ), hence 



* 7 (T) = - 



36C(3) 
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T-T. n 



(15) 



To verify our computational procedure, we also con- 
sider the ID case. The procedure that we used in 2D 
is also applicable to the ID case, with the only modi- 
fication that angular averages of the interaction in the 
particle- hole channel [Eq.®] are replaced by a sum of 
just two terms, for 9 = and 9 = tt, as the Fermi surface 
in ID consists of just two points k = ikp. The integrand 
in ([6]) vanishes at 9 = tt and diverges logarithmically at 
9 = 0, i.e., for D = 1, the renormalization of the interac- 
tion in the particle-hole channel also leads to logarithmic 



/ c (tt) = 2u - Uir, / s (7r) 



1 + 2uv In ^f 



The logarithmic corrections are cancelled out in / c (7r), 
but are present in J s (tt). The interplay between the 
behavior of the renormalized spin amplitude and the 
specific heat in ID is a more subtle issue, because the 
backscattering part of Sj(T) in ID contains extra 0(T) 
terms and does not reduce to Eq. ([T]) with the renormal- 
ized / s (7r) 



15, 24, 2 



To elucidate the difference between 2D and ID further, 
we consider a 2D system with a non-circular Fermi sur- 
face. Near an arbitrary point k.p on such a surface, the 
fermionic dispersion can be expanded as 



ek = v F k 



2m c 



(17) 



where fen and k± are the projections of vector k — k F on 
the normal and tangent to the Fermis surface at point k F , 
correspondingly, vf is the local value of the Fermi veloc- 
ity, and m c = 1/kvf is related to the local curvature, n, 
of the Fermi surface For a circular Fermi surface, 

k = kp so that m c = m = kp/vF and our 2D results 
are valid. If m c 3> m, the dispersion near two symmet- 
ric points ±k F is almost one-dimensional and, if only uq 
and are relevant, we should reproduce ID results. In- 
deed, evaluating the particle-hole contributions, labelled 
above as ((■■•)), for the case of to ^> m c , we find that 
they have the same logarithmic behavior, as in Eq. (fT())) : 
the only difference being that the logarithm is now cut 
by the largest of the two energies: T and E c = k F /2m c . 
Substituting these results into the expression for / c (7r), 
we find that to logarithmic accuracy 

J c (tt) = 2u -Un- 2(ttg + u\ - italic) \nmax(E c /T, 1) 

(18) 

For to c = 00, E c = 0, and the logarithmic term vanishes, 
just as it happens in ID. 

We see that the curvature of the Fermi surface is 



the crucial element of 2D consideration [2jJ. When the 
curvature is finite, the particlc-holc rcnormalizations of 
the scattering ampludes and of 6j(T) are not logatith- 
mic at the smallest T, and the logarithms only come 
from particle-particle renormalizations. Then f s (rr), 
f s (iT), and 8^f{T) are all logarithmically reduced. When 
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the curvature is zero, the dispersion is one-dimensional, 
particle-hole renormalizations also become logarithmic, 
and for / c (7r) the logarithms from the particle-particle 
and particle- hole renormalizations are cancelled out. In 
the eikonal approximation used by AE and in earlier 2D 
bosonization theories [28| , the curvature of the Fermi sur- 
face is neglected. In this situation, the charge amplitude 
behaves as in ID and is not renormalized, and S"f(T) re- 
mains linear in T. 

To summarize, in this paper we analyzed the effect of 
the Cooper-channel renormalization on the temperature 
dependence of the specific heat. We have shown the non- 
analytic term in the specific heat coefficient of a 2D Fermi 
liquid is expressed via the square of the full backscat- 
tering amplitude /(7r), renormalized in both particle- 
hole and particle-particle channels. Due to the particle- 
particle renormalization, both the charge and spin com- 
ponents of /(7r) are reduced by a factor of ln(Ep /T) in 
the limit of T —> 0. Consequently, the temperature de- 
pendence of 5j(T) is TS{T)/ln{E F /T), where S (T) 
is a slowly varying function, whose form depends on the 
details of the interaction in the particle-particle channel. 
When applied to ID, our method reproduces the can- 
cellation between particle-hole and particle-particle con- 
tributions to the charge channel. The logarithmic renor- 
malization of the charge amplitude is due to a specifically 
higher-dimensional effect-a finite curvature of the Fermi 
surface. For a flat Fermi surface, this renormalization is 
absent. 
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